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1 Introduction 



To begin with we are going to consider the inverse spectral problem for one-dimensional 
Schrodinger operator with periodic potential. In the late 60-s the famous discovery of the 
Inverse Scattering Transform for the KdV equation was done. A periodic analog of this 
transform was found in 1974. It is based on the solution of the following inverse spectral 
problem: 

to describe effectively the "isospectral manifold" of all the potentials with a given spectrum 
on the line (i.e. the spectrum of Schrodinger operator acting in the Hilbert space of square 
integrable complex-valued functions on the line R ). 

As everybody working in the quantum solid state physics knows, this spectrum gener- 
ically is a union of infinite number of intervals (allowed bands) on the energy line e. The 
complementary part on the energy line is also a union of infinite number of intervals (gaps 
or forbidden bands) whose lengths tend to zero for e — > +00. 

The periodic problem was solved in 1974-75 for the so-called finite-gap potentials that 
have only a finite number of gaps. Any periodic potential can be approximated by the 
finite-gap ones. This solution involves combination of the theory of Riemann surfaces and 
their ^-functions, Hamiltonian dynamics of special completely integrable systems and the 
spectral theory of Schrodinger operator. The mathematical technique used was (and still is) 
unusual for the community of physicists. Later the necessity to use this kind of mathematics 
appeared also in other branches of mathematical and theoretical physics (for example, in the 
string theory, matrix models, and supersymmetric Yang-Mills theory [3|, |J |5], [!], ^ 
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It seems to the authors this technique in future will be needed to 



the broad community of theoretical physicists. 

Integrability of the famous KdV equation u t = 6uu x + u xxx was discovered in 1965-68 (see 
ill) for the rapidly decreasing initial data on the line x. Exact solutions for the KdV 



equation expressing u(x,t) through the inverse scattering data of the Schrodinger operator 
L = —d^. + u(x, 0) were found. This procedure has been called the Inverse Scattering 
Transform (1ST). It was extended later for some other highly nontrivial (1 + l)-systems 
including such famous systems as Nonlinear Schrodinger NS± : iipt = —ipxx ± IV'PV' an d 
Sine(Sinh)-Gordon equations SG : u xt = sinu or u xt = sinhu. Note, that for the SG equation 
a large family of exact solutions was already constructed in XIX century by Bianchi, Lie and 
Backlund (see in |T9|, ^TJ). Beginning with 1974, several (2 + l)-dimensional physically 
interesting systems have been discovered as integrable by the 1ST procedure. The most 
famous of them is the KP system (see [p2] , f2~3|l). 

It is necessary to emphasize that the 1ST procedure in its original form can not be 
applied to the solution of the periodic problem (i.e. u(x,t) is periodic in the variable x). 
This problem was solved on the base of the new approach proposed in 
|25|, M, E7L H, |2|, |0, [31 (see the surveys 



[2 + l)-systems was found in pi R3, R5 



two-dimensional Schrodinger operator was started in 1976 (see 



24], in the works 
An extension of this method to 
New development of this approach associated with 

bsl m m m). 



HI) 



Complete detailed description of the solution of the periodic problem can be found in the 
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surveys, Encyclopedia articles ||40|| , and in the book [ 44 ] . We are going to present here basic 



ideas of this theory in the simplest form possible. Let us point out that KdV as well as other 
nontrivial completely integrable by 1ST PDE systems are indeed completely integrable in any 
reasonable sense for rapidly decreasing or periodic (quasi-periodic) boundary conditions, 
only. In fact, even that is well established for few of them. For example, for the KdV any 
periodic solution can be approximated by the finite-gap solutions. This statement easily 
follows from the theory of finite-gap potentials if we do not try to preserve the period, i.e. 
in the class of all quasi-periodic finite-gap potentials. The approximation of any periodic 
potential by the finite-gap potentials with exactly the same period, was constructed on the 
base of other approach developed in f46 |. The extension of the theory of Riemann surfaces 



and ^-functions to the specific class of surfaces of the infinite genus associated with periodic 



Schrodinger operator, was done in []4"? |. This theory is a beautiful description of the infinite 
limit. However, it seems that all fundamental properties of ^-functions associated with the 
complex continuation of variables, are lost in this limit. It is interesting to point out that 
analogous (but more complicated) theory of Riemann surfaces of the infinite genus was 
developed later in JO] for the periodic 2D Schrodinger operators. 

Outside these functional classes almost no effective information is known. Beautiful 
methods have been developed also for the studies of the special self-similar and "string- 
type" solutions, but in most cases they lead to the very hard analytical problems associated 
with the famous Painleve equations and their generalizations (E8L |49|, [50, [y], E3l 1521 154], |55|). 



2 Rapidly decreasing potentials and GGKM Proce- 
dure. Backhand Transformations 

Let us recall the basic information about the 1ST method for KdV. We start from the so-called 
Lax Representation for KdV (see The Heisenberg type equation for the Schrodinger 

operator L 

L t = [L, A] = LA - AL, (2.1) 
L = -d 2 x + u, A — — 49| + 3(ud x + d x u), (2.2) 
is equivalent to the identity (KdV equation) 

u t = Quu x + u xxx . (2.3) 

For this reason, any KdV type equations admitting some analog of the Lax representation is 
called isospectral deformations. The existence of such deformations indicates the possibility 
of effective solution of the Inverse Scattering Problem for the operator L. 

For the rapidly enough decreasing functions u(x, t) — > 0, x — > ±oo we define two bases of 
solutions (t is fixed): 

(f)±(x,t; k) ~ exp ±lkx , x — > -oo, L<fi± = A0±, (2.4) 
ip±(x,t; k) ~ exp ±ikx , x -> +oo, L?p± = \ip±, k 2 = A. (2.5) 
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By definition, monodromy matrix T connects these two bases, Tcf) = tjj for the column vectors 

(f) = ip = 



T = , , ip+ = a( t ) + + b(j)-, ip- = c0 + + . (2.6) 




A conservation of the Wronskian implies that detT = ad — be = 1. For the real values of 
k or A > we have a = d,c = b. Therefore, \a\ 2 — \b\ 2 = 1. The whole set of the so-called 
inverse scattering data can be extracted from the monodromy matrix T if it is well-defined 
for all complex values of k. The so-called scattering matrix is constructed from T for the 
real k. Its entries are the transmission coefficient 1/a and the reflection coefficient b/a. The 
property 6 = for all real values of k characterizes reflexionless or multisoliton potentials. 
For all rapidly decreasing potentials the matrix element a(k) is well-defined for complex k 
such that Im k > and a — > 1 for k — > oo, Im k > 0. There is only a finite number of 
purely imaginary zeroes a(k n ) = in this domain. They correspond to the discrete spectrum 
A n = kl < 0. 
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The famous result of |17| (GGKM procedure) easily follows from the Lax representation, 
which implies the equation: 

T ' = I T - A )' A =("o* 4*) (") 

This result was formulated as a set of the following GGKM formulas 

at = 0, bt = -8(ik) 3 = -ct, d t = 0. (2.8) 

The latter equations give a full description of KdV dynamics in these variables because any 
rapidly decreasing potential can be reconstructed from the inverse scattering data. A special 
family of the reflexionless potentials where 6 = for real values of k, leads to the so-called 
multisoliton solutions for KdV equation (see [13].) 



The multisoliton solutions can be also directly obtained with the help of the elementary 
substitutions (Baclund Transformations) transforming any solution of the KdV into another 
solution: let u be a solution of the KdV equation and v be a solution of the Ricatti equation 
a + u = v x + v 2 with the initial value independent of time. The new function u = —v x + v 2 
satisfies the KdV equation. Starting from the trivial solution u = uq = we construct 
a sequence of potentials u n -i = u n ,n > 0, given by the Biicklund Transformation. We 
choose parameters a n , a\ > a<i > ... > a n > . . ., and take the real nonzero functions 
/„ -> oo, x -> ±oo, -f nxx + u n -xfn = oc n f n , which define v n = (}ogf n ) x . Every such 
sequence leads to the multisoliton reflexionless potential: 

2a 

Uo = °' Ui = ~ ch^(x - *o) + wr • • • (2 - 9) 

In terms of the Schrodinger operator this transformation (invented by L.Euler in 1742) is 
called Darboux Transformation. The operator L can be factorized 

L = -d 2 + u = -(d + v)(d - v). 



4 



Using the non-commutativity of these factors, we define the Darboux Transformation for the 
operator and its eigenfunction in the following way: 

L = —d 2 + u = —(d — v)(d + v), i> = (d + v)ip . (2.10) 

These transformations can be considered as some kind of discrete spectral symmetries for 
Schrodinger operators. They preserve a spectrum of the operator L (maybe except for one 
eigenfunction). 



3 KdV Hierarchy. Integrals of motion. Hamiltonian 
formalism. 

The local integrals for the KdV equation can be constructed with the help of the Schrodinger 
operator. Consider the associated Riccati equation v x + v 2 = u — k 2 and find the solution 
for it as a formal series in the variable k: 



v[x 



,k) = tk + j2v n (x)(ik)- n , (3.1; 



n=l 



where all v n are polynomials in the variables u, u x , . . . The integral along the line x is a 
fc-dependent constant of motion for the KdV equation 

For the real potentials u(x,t) and real k we can see that the imaginary part of v(x, k) — ik 
is a total derivative. The remaining quantities in the expansion : 

J (v(x, k) — ik)dx = ^2 J v n (x)(ik)~ n dx., (3.3) 



n>l 



define local integrals of motion 

I n = c n J v 2n +3(x)dx, n = -1, 0, 1, 2, . . . , (3.4) 

where constants. After the proper choice of the constants c n we have 

I-i = Judx,I = Ju 2 dx,I\ = J(u x /2 + u 3 )dx, . . . , (3.5) 

I nt = 0, u t = Quu x + u xxx . (3.6) 
Let us introduce a GZF Poisson bracket |56|, |57|] on the space of functions 

{u(x),u(y)} = 5'(x — y)). (3.7) 

Then any functional H (Hamiltonian) defines the corresponding Hamiltonian system 



5 



For the case H = J_i we get a trivial flow (i.e. this integral is a Casimir for the GZF 
bracket). For H = I Q we are coming to the ^-translations u t = u x . Let us call this equation 
KdVo. For the case H = Ii we have the ordinary KdV = KdV\. Higher integrals give us 
the equations KdV n of the order 2n + 1 admitting the Lax representations with the same 
Schrodinger operator L but with the differential operators A n = (const)d x Vn+1 + . . . : 

^ =a -{m) =v - u (3J) 

In particular, A = d x , A\ — A. Nice formula for all operators A n can be extracted from 



Let L = —C 2 , where C = d x + Y^k>i a k(u, u x , ■ ■ -)d~ k . Here all are polynomials in the 
variables u, u x , . . . and d~ l a = J2 n >o(~ ^-) n 0'^d~ n_1 for the composition of the operator d^ 1 
and multiplication operator by a. By definition, 

A n = (C 2n+1 ) + = (L 2 " +1 / 2 ) + , (3.10) 

where the sign + means omitting of all strictly negative powers of d x . 

All higher KdV n systems can be integrated by the same 1ST procedure for the class of 
rapidly decreasing functions. In particular, GGKM equations for the scattering data (or 
monodromy matrix) have the form 

T tn = [T, A n ] , A n = (canst) ( { _^ )2n+1 ) - (3.11) 

The latter result also implies that all these flows commute with each other. Hence, we get 
the following conclusion without any calculation: 

integrals I n have zero Poisson brackets, 

J ou[x) ou(x) 

A generalization of the GZF Poisson bracket for the isospectral deformations of the higher 
order (scalar) Lax operators L was found in |59| . 

It should be emphasized that there exists a family of local field-theoretical Poisson brack- 



ets (LM-brackets, |60[) describing the KdV theory: 

B = B M = Xd x + n(-d% + Aud x + 2u x ) = \B + ^B u (3.13) 

{u(x),u(y)}x lfl = B8(x-y), (3.14) 

{4,AnW = 0- (3.15) 



These brackets were generalized for the higher order operators L in |pl . 

The recurrence operator B 1 B$ 1 = C generates all the right-hand sides of all higher KdV 
systems: 

C(0) = u x , C 2 (0) = 6uu x - u xxx , C n (0) = d x f^yl • (3-16) 
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It gives also a simple proof of a very useful identity (|62|): 
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dx = (const) I n -\. (3-17) 



5u(x 



All these identities are local and can be used for the class of periodic functions as well. 
However, we shall see in the next section that the direct analog of GGKM procedure does 
not lead to the integration procedure. We are going to use a different approach. 



4 Spectral theory of periodic Schrodinger operators. 
Finite-gap potentials 

The spectral theory for the periodic potentials on the whole line x is based on the monodromy 
matrix as in case of the scattering theory. However, in the periodic case with a period 
T < oo, we have nothing like the selected point x = oo for the definition of the monodromy 
matrix (as it was for the rapidly decreasing case (T = oo)-see Section 2 above). Any 
point xo can be used. Let us fix initial point x$ and choose a special basis of the solutions 
C(x, xo, e), S(x, xo, e) for the spectral equation LC = eC, LS = eS such that for x = xq we 
have: 



C S 




(4.1) 



The shift operator T : x — > x + T in the basis C, S defines the monodromy matrix 

?(xo,e)=(l b A, (4.2) 



C(x + T) = a C{x) + b S{x), S(x + T)=c C(x) + d S{x). (4.3) 

The key element of the periodic spectral theory is a notion of the, so-called, Bloch waves or 
Bloch-Floquet eigenfunctions. We present here some essential properties of these functions 
without any proofs. An exposition of this theory may be found in the Encyclopedia article 
|40fl , where the main ideas of the proofs are clearly presented for the difference Schrodinger 



operator (it is much simplier). 

By definition, the Bloch-Floquet functions are solutions of the Schrodinger equation that 
are at the same time eigenfunctions of the shift operator, i.e. 

Lip = eip, Tip(x) = ip(x + T,e) = exp(±ip(e)T)ip (x) . (4.4) 

We uniquely normalize ip by the condition ip\ x=xo = 1. 

For any complex number e the eigenvalues, w±(e) = exp(±p(e)T), of the shift operator 
are defined by the characteristic equation for the monodromy matrix. ^From the Wronskian 
property it follows that det T = 1. Therefore, the characteristic equation has the form 

w 2 - (tr f)w + 1 = 0. (4.5) 
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Figure 1: 



The multivalued function p(e) is called quasi-momentum. 

The spectrum of the Schrodinger operator on the whole line is a union of spectral zones 
which are segments of the real line of the variable e, where the quasi-momentum is real. The 
latter condition is equivalent to the inequality: 

|Tr T\ = 2cos(pT) < 2. (4.6) 

The typical graph of the function /(e) = cos(pT) is presented on the fig 1. In particular, its 
extreme points f'(e) = "generically" are located inside of the gaps (i.e. for the open and 
everywhere dense set of periodic potentials we have |/| > 1 at the extreme points, and there 
is only one extremal point in each gap). For some special cases we may have / = ±1 at 
the extreme point. Such point lies inside the spectral zone. However, generic perturbations 
create new small gap nearby this point — see fig 2. This point is a double point of periodic 
or anti-periodic spectrum with the boundary conditions i]){x) = ±tp(x + T). 
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Figure 2: 
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The Riemann surface of the Bloch-Floquet functions is defined by the equation 

z 2 = cos(p(e)T), (4.7) 

but this surface is nonsingular only for the generic case (when there are no double points of 
the periodic or anti-periodic problem for the Schrodinger operator). For the large e — > +00 
the following asymptotics is valid for the gaps (see fig. 3). 

1. The length of gaps tends to zero with its rate depending on the smoothness of the 
potential (this rate is exponential for the analytic potentials). 

2. All gaps are located nearby the points e m = 4tt 2 m 2 T~ 2 , and the distance of gap from 
this point tends to zero. 

A nonsingular Riemann surface of the Bloch-Floquet solutions is defined with the help 
of the equation 

z 2 = (e-e )(e-e 1 )...(e-e n )..., (4.8) 

where eo < e\ < . . . , are simple eigenvalues of the periodic and anti-periodic spectral prob- 
lems. 
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Figure 3: For the generic T-periodic potentials the gaps are located near the 
points e m = 47im 2 T~ 2 for e — > 00. Here m are the integer numbers. 



As it follows from the Lax representation for all higher KdV n systems (|3.9| ), these bound- 
ary spectral points for the periodic operator L in the Hilbert space L 2 (R) of the square 
integrable complex- valued functions on x-line are integrals of motion for the KdV hierarchy. 

We can say that the Riemann surface of the Bloch-Floquet solutions as a whole is an 
integral of motion for the KdV hierarchy. How many potentials correspond to the same 
Riemann surface (i.e. have the same spectrum in L 2 (R))1 

The original idea to introduce a special class of potentials for which this problem can 
be solved effectively was proposed in p4[ and is based on the KdV hierarchy. This idea 



naturally combines two different ways to describe the corresponding isospectral manifold of 
potentials. Both are in fact closely related to each other. 
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The first approach — use of the KdV hierarchy and Hamiltonian dynamics: 

Let us consider a stationary equation for some linear combination of higher KdV n flows 
( p.9| ). It is an ordinary differential equation that can be written in the form: 

5(I n + ciVi + . . . + c n I + c n+1 J_i) = 0. (4.9) 

This is a finite-dimensional Hamiltonian system with n degrees of freedom depending on 
n + 1 parameters ci, . . . ,c n+ \. It is a completely integrable Hamiltonian system because 
there is a pencil of the commuting flows. This pencil coincides with the restriction of all 
higher KdV m on this stationary subset of functions given by the equation ( f4.9| ). Therefore, 
its generic nonsingular solution is expected to be a periodic or quasiperiodic function of x. 

In the next section we shall construct for (|4.9| ) some kind of the Lax representation 

n+1 

A* = [Q(x, A), A(x, A)], A = A„, + ]T c fc A fc , (4.10) 

k=l 

using 2x2 traceless matrices depending on the parameter A and (u, u x , . . .) polynomially. In 
the today terminology some people call them loop groups. Note that in all soliton systems this 
A-loops have very specific A-dependence (polynomial, rational and in some exotic examples 
- elliptic functions) . 

Lax representation ( |4.10|) implies that an algebraic Riemann surface defined by charac- 
teristic equation 

det(A(A) - zl) = P(A, z) = 0, (4.11) 

does not depend on x and therefore is an integral of ( |4.9|) . The same Riemann surface can be 
extracted from the commutativity equation [L,A] = 0, which according to ( |3.9| ) is equivalent 
to (0). 

It should be emphasized that the latter form of (|4.9| ), i.e. the commutativity condition 
for two ordinary differential operators 

[L, A n + Y, c k A n . k ] = [L, A] = 0, L = -d 2 x + u, (4.12) 

was considered formally (i.e. locally in the variable x without any periodicity assumptions) 
as a pure algebraic problem in 1920-s (see |53], 0). Even the formal algebraic Riemann 
surface (|4.11|) appeared as a relation P(L, A) = 0. 

According to our logic however, the corresponding system of equations in the variable x is 
Hamiltonian and completely integrable. Therefore, its generic solution is quasiperiodic in x, 
containing a dense family of periodic solutions. We may ask about the spectrum of the cor- 
responding operators L in L 2 (R) and boundaries of gaps. Remarkably, they exactly coincide 
with the branching points of the Riemann surface defined by equation (|4.11| ). Therefore, for 
the periodic potential which satisfy equation ( |4.9| ) the nonsingular spectral Riemann surface 
of the Bloch-Floquet solutions is an algebraic Riemann surface of genus g = n. The spectrum 
of such an operator contains only a finite number of gaps [e^j-x, €23}, j — 1, • • • ,n. 

This key step unifies the first approach with the second one (below). So the solution of 
the inverse spectral problem can be identified with the process of solution of some special 
families of completely integrable systems using Riemann surfaces. 
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Figure 4: The endpoints of gaps 6j are the branching points of the 2-sheeted 
Riemann surface. The poles jj are located in the gaps on one branch only. (Here 
n = 3) 



The second approach — periodic spectral theory in the Hilbert space L 2 (R): 

For any periodic Schrodinger operator L = —d 2 + u(x = t ,ti,...) we have already 
defined the Riemann surface T of Bloch-Floquet solutions with the help of the monodromy 
matrix T. For the finite-gap potentials a graph of the function / = cos(p(e)T = l/2tr T is 
highly degenerate (see fig. 2 and compare with fig.l). For all real and large enough e we have 

I/I<1- 

Once again, we ask how to describe all the potentials with the same finite-gap spectrum 
and what are additional variables which uniquely define a finite-gap potential? It turns out 
that these additional spectral data are the poles of the Bloch-Floquet functions. It can by 
shown that the Bloch-Floquet solutions defined by Q4.4Q and normalized by the condition 
ip\ x =x = 1 have exactly one simple pole jj(xo) inside each gap or on its boundary. A 
point of the hyperelliptic surface T ( |4.8|) can be represented as the complex number e and 
a sign (or a branch) of the radical z = JJe — eo)(e — t\) • • •■ The branches of the radical 
coincide at the boundaries of gaps. Therefore, each gap should be considered as a closed 
cycle cij,j = 1, 2, ... n (see fig. 4). Any pole of ip should be considered as a point jj on the 
cycle dj. The above-mentioned statement that there is only one pole in a gap means, in 
particular, that ip h as no pole at the point (7^, +) if it already has pole at (7^, — ) and vice 
versa. Geometrically, the total set of poles (71, . . . ,7 n ) represents a point on the real torus 
T 9 = ax x . . . x a n . As we shall see later, this set of poles completely determines the original 
potential. Different points of the real torus lead to the different potentials (if normalization 
point x is fixed). 



5 Periodic analog of GGKM. Zero Curvature repre- 
sentation for the KdV hierarchy and corollaries 

As it was emphasized above, the definition of the monodromy matrix in the periodic case 
depends on a choice of initial point Xq. For different choices of the initial point the cor- 
responding monodromy matrices are conjugated (because they represent the same linear 
transformation in different bases). Therefore, a dependence of the monodromy matrix with 
respect to the choice of the initial point xq can be described by the equation: 

r„ = [QM„),r], Q=( £ _° w J). (6.1) 
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In the same way, for the isospectral deformations corresponding to the KdV hierarchy (i.e, 
all higher KdV n systems) we can establish the following equations (the periodic analogs of 
GGKM): 

T tj = [Aj, T], A = Q, t = x . (5.2) 

A compatibility condition of equations ( |5.2j ) for any pair of variables ti, tj implies the following 
zero curvature representation for the KdV hierarchy (where periodic boundary conditions are 
already inessential) 

[d ti -A i ,d tj -A j } = 0. (5.3) 

Using the Lax representations ( p.9|) for all higher KdV n , we can express all the matrices Aj as 
polynomials in the variable e and variables u(xq),u xo (xq), . . .. For example, for the ordinary 
KdV we have (replacing x by x): 

A - I ~ Ux 2w + 4e \ , , 

Al ~ I -4e 2 + 2tu + 2u 2 - u xx u x J 



The matrices A k can be completely reconstructed from equations ( |5.3|) and from the following 
properties of matrix elements: 

A k = (const) \ ak bk ) , (5.5) 
\ Ck d k ) 

d k + a k = TrA fc = 0, b k = e k + u^e^ 1 + ..., (5.6) 

det A k = -(a 2 + be) = (const) R 2k+ i(e) = (const)(e 2k+1 + ...), (5.7) 

2ctfc = — b kx , (i?2fc+i)x = —b k - (5.8) 

For i = 0, equations ( |5.3| ), where t = x, give the zero-curvature representation for the KdV n 
system: 

d x A n -d tn Q = [Q,A n )- (5.9) 

From this representation and from the periodic analog of GGKM we are coming to the 
following results for the stationary higher KdV n system ( (4.10| ): 

1. The monodromy matrix T commutes with A 

[T,A] = 0. (5.10) 

Therefore, they have common eigenvectors. It implies, in particular, that the Bloch-Floquet 
function ifi± is single-valued on the algebraic Riemann surface associated with the matrix A; 

2. The stationary higher KdV admits an e-parametric Lax-type representation in the 
variable x ( ^.10| ): 

(A) X = [Q,A], A = A n + Y / c k A k . (5.11) 

Therefore, we have a full set of conservation laws organized in form of the Riemann surface 
(i.e. all coefficients of the polynomial P are re-independent): 

det(A(e) - zl) = P(e, z) = 0, (5.12) 



12 



P(e,z)=z -R 2n+ i(e) = (const) (e - e ) ...(e-e 2n )- (5.13) 

These points €j are exactly the boundaries of gaps for periodic potentials because the Rie- 
mann surfaces of the matrices T and A coincide. The matrix element &12 for the matrix A 
determines another set of points 

612(e) = {const){t - 7l ) . . . (e - 7„). (5.14) 

These points coincide with projections of the zeroes of Bloch wave ip± as the functions of 
the variable x or poles as the functions of Xq (see the next section). 

An original approach to the solution of the inverse spectral problem was based on the 
use of the following trace type formula for the potential: 

— u(x)/2 + const = 7i(x) + . . . j n (x). (5.15) 

iFiom ( |5.11| ) and ( |5.14| ), "Dubrovin equations" defining dynamics in x of the points jj can 
be derived. They have the form 



7ix = v/7ra7j \ (5.16) 

These equations (see |27|, [32]) can be linearized by the, so-called, Abel transformation. 

Consider the first kind differentials on the Riemann surface (i.e. holomorphic 1-forms 
without poles anywhere - even at the infinity). The basic first kind forms are 

e j ~ 1 de 

Vj = , j = 1, ■ ■ ■ ,n. (5.17) 

lR 2 n+i(e) 



It it is convenient to choose a normalized basis taking linear combinations Q s = J2j v sj^j 
such that 

fl, = 5e 



S]1 

■J 



for all the gaps cij. 

Fix a set of paths Kj from the point Pq = 00 to the points Pj on the Riemann surface. 
Abel transformation is defined by the formula 

A p (P u ...,P n )=J2 I % ■ (5-18) 
Equations ( ^.16|) after the Abel transform become linear: 

A px = U q = ±- I dp . (5.19) 

Z7T Jbg 

Here the closed paths b q are "canonically conjugated" to the paths aj (see fig. 5). It means 
that the intersection numbers are aj ob q = 5j q . The differential dp is equal to p E de where p(e) 
is a multi-valued quasimomentum. It is a second kind differential (i.e. meromorphic 1-form 
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Figure 5: The canonical basis of cycles {ai,bj) on the Riemann surface T with 
branching points e , . . . , e 2n - (Here n — 3) 



on the surface T with a pole of order 2 at the point oo with negative part dk = —(dw)w 2 
in the local coordinate w = k~ l near oo), and therefore has the form 

fr^ + W 1 - 1 -')*, (5 . 20) 
yi?2n+i(e) 

All the coefficients Vj can be found from the normalization condition: 

' dp = 0, j = l,...,n. (5.21) 

Following the classical XIX century theory of Riemann surfaces, formulas ( |5.15| - |5T21"|) lead 
to some expression of the potential through the ^-functions avoiding the calculation of the 
eigenfunction ip± (see [TZI| [27], |29|, [JTJ). The most beautiful functional formula for the 
potential was obtained in the work [^9] . 

However, in the next section we shall not follow this approach for solving the inverse 
problem which has been outlined above, but use another approach proposed for the first 
time in the article |3^] (see appendix, based on the idea of A. Its): it is possible to calculate 
all the family of the Bloch-Floquet functions. The original approach has been used later in 
the problems, where the effective calculation of eigenfunction is impossible (like in the cases 



of higher rank commuting OD linear operators (see for example [37]). 



6 Solution of the periodic and quasiperiodic inverse 
spectral problems. Baker- Akhiezer functions 

We are going to solve the inverse spectral problem for the finite-gap potentials following the 
scheme proposed in [34], |35| and based on the concept of the Baker-Akhiezer functions. 
These functions are uniquely defined by their analytical properties on the spectral Riemann 
surface. As we shall see later this scheme is evenly applicable to the solution of inverse 
problems in two-dimensional case where the corresponding analytical properties naturely 
generalize the analytical properties of the Bloch-Floquet solutions for finite-gap Schrodinger 
operators. 

Let us start with the following real inverse spectral data: 
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1. Riemann surface T given in the form 



z 2 = (e - e ) . . . (e - e 2n ), 

where all numbers ej are real; 

2. a set of n real points 7j G T such that there is exactly one point jj on the cycle a.,-. 

Below we shall for brevity identify the point 7j as a point on T with its projection on 
complex e-plane. The condition that there is only one point on each cycle means that 
projection of the points satisfy the restriction: 

e 2 j-i < lj < e 2j , j = 1, ... 71, (6.1) 

By definition, the complex inverse spectral data are the same data where the Riemann surface 
(i.e. its branching points e^) and the points jj G T are arbitrary complex points. 

As we shall see below a generic set of algebraic-geometric spectral data leads to the 
explicit solution of the inverse spectral problem in terms of the Riemann t/ieta-function. The 
corresponding potentials are complex meromorphic quasi-periodic functions of the variable x. 
For the real data described above, we are coming to the smooth (even analytic) quasiperiodic 
potentials with their periods expressed through the hyperelliptic integrals (see below). It is 
necessary to mention that there is no way to find simple and effective criteria for the potential 
to be periodic in terms of these data. The periods depend on the Riemann surface, only. 
Of course, we may write the condition that all the corresponding hyperelliptic integrals are 
commensurable, however this condition in transcendental. Recently, based on the results of 
|7], ^5] , the effective numerical approach has been developed for the solution of this problem. 



It is based on the discovery of some specific dynamic systems on the set of potentials which 



preserve all the periods but change spectrum. Following []65[ ) we start from one periodic 
potential and create all others using these dynamic systems. 

In the last section we shall present some examples of finite-gap potentials written in terms 
of elliptic functions. They are periodic in the variable x (even double-periodic as functions 
of the complex variable x). The first nontrivial examples different from the classical Lame 
potentials u(x) = n{n + l)p(x) were found in fl25"f . This subject was developed in the works 

We define the Baker- Akhiezer function ip = ip(x, ti, . . . t n ; P) for the parameters x, tj and 
the point P = (e, ±) G T, by its analytical properties on T with respect to the variable P. 
For the case ti — ty, — ■ ■ ■ — these analytical properties are just the same as the analytical 
properties of the Bloch-Floquet solutions of the periodic finite-gap operator. 

^From pure algebraic-geometric arguments it follows that there exists a unique function 
■0 such that 

1. it is meromorphic on T outside the infinity and has at most simple poles at the points 
lj, j = l,...,n; 

2. in a neighborhood of the infinity the function if) has the form 

i/j = exp[xk + t x k 3 + . . . + k 2n+1 t n }(l + £ x (x, t)^ 1 + ...) , 
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where k 2 = e and therefore, k 1 is a local coordinate on the Riemann surface near infinity. 

A proof of this statement is identical to the proof of the existence and uniqueness of gen- 
eral Baker- Akhiezer functions that were introduced in for construction of exact solutions 



of two-dimensional KP equation and all associated Zakharov-Shabat hierarchies. 

General Baker- Akhiezer function is defined with the help of an arbitrary Riemann surface 
of finite genus n instead of special hyperelliptic (i.e. two-sheeted) surfaces. We fix an 
arbitrary "infinity" point Pq on it; local coordinate k~ x = w where w(Pq) = ; generic set of 
points (71, . . . , 7„), and numbers (x = n, r 2 , . . . , rfc, . . .). The corresponding Baker-Akhiezer 
function has the same analytical properties as above but, (2) is replaced by the following: 

V> = exp[xA; + nk 2 + r 2 k 3 + . . .](1 + ^(x, r)^ 1 + ...). (6.2) 

For the KdV hierarchy we have a hyperelliptic Riemann surface, k 2 = e, and r 2 k = 0, T 2 j + i = 
t r 

The use of such type of functions (dependent only on a single parameter x) has been 



proposed by Baker in his note ||71| , for the common eigenfunctions of two commuting OD 
linear operators. He expected that this construction will improve the results of [^, |6^| . He 
made also very interesting conjecture that this approach may seriously improve the classical 
theory of ^-functions. Unfortunately, this program was not realized and was forgotten. The 
soliton theory appeared many decades later absolutely independently. In the 70-s it started 
to use such kind of functions on Riemann surfaces in the process of solution of a periodic 
problem for the KdV type systems and inverse spectral periodic problems. In the classical 



spectral theory Akhiezer (|[72||) was the first to use some special cases of this function for the 
construction of some examples of operators on the half-line x > with interesting spectral 
properties. No one of the authors before the 70-s had associated anything like that with 
periodic problems. Indeed, this type of functional construction on the Riemann surface was 
extracted in 1974 from the work 1721 . 



We shall prove the existence and uniqueness of the Baker-Akhiezer function and present 
its exact expression through the ^-functions later. At this moment we would like to show 
how the uniqueness of this function leads to the proof that ip is an eigenfunction for the 
Schrodinger operator for our special case. Apply twice the operator d = d x to this function 
and use analytical properties. After elementary calculation, we are coming to the formulas: 

dip = kip + exp[kx + . . .\iZixk~ 1 + •••), 

cPtP = k 2 tp + 2£ lx ip + 0(k- 1 )exp{kx + ...]. (6.3) 
^From the latter equations, we get the equality 

(d 2 - k 2 - 2^)^ = Oik' 1 ) exp[kx + ...]. 

The left-hand side is globally well-defined function on the same Riemann surface because 
k 2 = e. It has the same poles independent of parameters. Up to the same exponential factor, 
it is of the order 0(fc _1 ) at the infinity. Therefore, it is equal to zero due to the uniqueness 
of the Baker-Akhiezer function. 
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So we are coming to the conclusion that 



Lip = eip, L — -d 2 + u, u = -2£ lx . (6.4) 

We can apply the operators dj = d tj to the function ip: 

djip = k 2j+1 ip + exp[kx + . . -K^ij/r 1 + &jk- 2 ). (6.5) 

As in the previous case, we can easily construct a linear operator Aj = d 2 ^ +1 + . . . with the 
coefficients independent of k such that (dj — A n )ip = Oik^ 1 ) exp(kx + ...). Using global 
analytical properties on the Riemann surface, we deduce from this that the left-hand side is 
equal to zero as before. The compatibility condition of these pairs of equations is exactly a 
.fTc/V^-system for the potential u(x, t\, . . . , t m ). 

Following the works |33|, [34], [35] , we can prove in the same way that for the general 
Baker- Akhiezer function ip(x, r, P) associated with arbitrary Riemann surface with fixed 
local coordinate near a puncture, the following equations are valid: 

(d Tk -L k )ip = 0, k = 2,3, (6.6) 

Here = d% + . . . are linear OD operators acting on the variable x with coefficients de- 
pending on parameters r. These coefficients are differential polynomials in the coefficients of 
the expansion of the regular factor of the Baker- Akhiezer function at the puncture (infinity). 
They are uniquely defined by the condition that the expansion of the left-hand side of ( |6.6j ) 
at the puncture has the form 0(k~ 1 )ip. 

The compatibility conditions 

[d Ti -L j ,d Th -L k ]=U. (6.7) 

are equivalent to non-linear partial differential equations for the coefficients of the operators 
Lfc. For the case j = 2, k = 3, r 2 = y, r 3 = z, the operators L 2 and L 3 have the form 

3 

L 2 = d 2 x - u(x, y, t), Lz = dl- -ud x + w(x, y, t), (6.8) 

where we consider the dependence of the coefficients with respect to the first three variables 
Ti = x, r 2 = y, r 3 = t, only . The coefficient u(x, y, t) is equal to u = 2£i x (x, y, t), where £i is 
the first coefficient of the expansion ( |6.2j ) . 

^From (|6~7j ) we get a system of two equations for two coefficients u(x, y, t) and w(x, y, t) 
which can be reduced for an equation for u(x, y, t). The reduced equation is the famous KP 
equation 

3u yy = (Au t - Quu x + u xxx ) x . (6.9) 

It appeared in physics literature for the investigation of the transversal stability of the KdV 
solitons (see | |73]] ) and is one of the most natural physical two-dimensional analogs of KdV. 



Lax representation for it was found in |22], We call the whole set of higher systems (|6.7| ) 
the KP hierarchy. We shall discuss the periodic problem for this equation at greater length 
in the next section. 
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For the special choice of the hyperelliptic Riemann surface and k = ±^/e we have k 2s = e s , 
where the function e is well-defined globally as a meromorphic function on the Riemann 
surface T. Therefore, we may represent globally the corresponding Baker-Akhiezer function 
in the form 

ifj = exp{r 2 k 2 + r A k 4 + ...), (6.10) 

where ip does not depend on the parameters r 2 j. So in this case all the KP hierarchy reduces 
to the KdV hierarchy. 

The relationship of this constructions with commuting OD linear operators is as follows. 
Let f(P) be a meromorphic function on the Riemann surface T with one pole at the puncture 
Pq. Its negative part written in the parameter k is some polynomial q(k) — q[ + q 2 k l_1 + 
. . . + q\k. Apply the operator Af — q\L\ + q 2 Li_i + . . . + q\d x to the function vp. By the 
definition of the operators L^, we can see that 

(Af - f)ip = 0(k~ 1 )ip. (6.11) 

We conclude as before, that the difference Afip — f(P)ip = 0. 

For any pair of functions / and g on the Riemann surface T with poles at the point oo, 
we get a pair of the commuting OD scalar linear operators Af,A g such that AfA g = A g Af. 
In a special case of the hyperelliptic Riemann surface z 2 = (e — eo) . . . (e — e 2n ), we have pair 
of functions f — e, g — z, leading to the Schrodinger operator L = —d 2 + u commuting with 
the second operator of order 2n + 1, because z = k 2n+1 + . . . 

Now we return back to the problem of existence and uniqueness of the Baker-Akhiezer 
function. The simplest way to prove this existence is to define this function by exact for- 
mula in terms of the ^-function and meromorphic differentials. Let us recall first necessary 
information. 

As in the previously discussed case of hyperelliptic curves, we introduce a basis of cycles 
a,, bj,j — 1, . . . , g, on a Riemann surface T of genus g with canonical matrix of intersections: 
aj • bj = 5ij, and a basis u>i of holomorphic differentials normalized by the condition 

j> ujj = Sij. (6-12) 

The matrix B = (B>ij) of 6-periods of these differentials 

Bij = <f> Uj, (6.13) 

J hi 

is symmetric and has positively defined imaginary part. The Riemann ^-function is a function 
defined with the help of this matrix by the formula: 

6{z\B) = e 2niiz ' m)+{Bm ' m) (6.14) 

where z — (z±, . . . , z g ) is a complex g- dimensional vector, (m, z) stands for a standard scalar 
product and summation is taken over all integer vectors m = (mi, . . . ,m g ). Theta- function 
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is an entire periodic function of g variables Zj and has the following monodromy properties 
with respect to the shifts defined by vectors B k which are columns of the matrix of 6-periods: 

9{z + B k ) = e - 2 * iz "-* iB '"'6(z). (6.15) 

The basic vectors e k and the vectors B k define a lattice C in C 9 which determines g- 
dimensional complex torus J{T) = C 9 /C called Jacobian of the curve. The Abel map 
A : T — > J(r) is defined by the formula 

'■i> 



A k (P) = f u k . (6.16) 



Note that the vector A(P) with coordinates A k (P) depends on the choice of path of integra- 
tion but its ambiguity just coincides with shifts by vectors of the lattice C 

iFiom. the monodromy properties of ^-function it follows that zeros of the multivalued 
function 6(A(P) + Z) considered as a function on T are well-defined. For a generic vector 
Z this function has exactly g zeros (71, . . . , 7 9 ). The vector Z can be expressed in terms of 
Abel transforms of these points by the formula 

3=1 

where K is a vector of the Riemann constants. 

Let us introduce a set of meromorphic differentials dfli that are holomorphic on T outside 
the puncture where they have poles of the form 

cKl i = dk i (l + 0(k- i - 1 )), (6.18) 

and normalized in a usual way by the condition 

dVLj = 0. (6.19) 

J at 

The Abelian integrals 

fi.(P) = J dVti (6.20) 

are multivalued functions on Y. 

Let Uj be a vector with the coordinates 

u »=-kL ia >- (6 - 2l) 

Then from the statements presented above, it follows directly that the formula 

p ) = "P(S^(P)) e {M p) + zmzMt, + z) ' (6 ' 22) 

correctly defines a function on T which satisfies all the properties of the Baker-Akhiezer 
function. 
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Suppose now that there exists another Baker- Akhiezer function ipi. ^From the definition 
the Baker- Akhiezer functions it follows that the ratio ipi/ip is a meromorphic function on T 
which is equal to 1 at the puncture and with only possible poles at the zeros of the function 
tp. According to (|6.22j ) the zeros of ip are zeros of the function 9(A(P) + J2iUiTi + Z). 



Therefore, ip has g zeros. The simplest form of the Riemann-Roch theorem (which can be 
considered as generalization of the Liouville theorem for Riemann surfaces) implies that a 
function on T with at most g poles at a generic set of points is a constant. Therefore, ipi = ifi 
and the existence and uniqueness of the Baker- Akhiezer function is proved. 

Now according to the previously established formula u = 2d x £\, in order to get an exact 
formula for a solution of the KP hierarchy it is enough to take the first coefficient £i(t) of 



the expansion of the pre-exponential factor in (|6.2j ) at the puncture. Finally we obtain the 
expression 

u(t) = 2d 2 x log 0(52 U i T i + z ) + const ( 6 - 23 ) 

i 

for the finite gap solutions of the whole KP hierarchy. 

If we consider only the KP equation, we get the formula 

u(x, y, t) = 2d 2 x log 6(Ux + Vy + Wt + Z) + const, (6.24) 

where we redenote x = Ti, y = r 2 , t = r 3 and U = U±, V = U2, W = U3. 

For the case of the hyperelliptic curve the vector V = and we get the Its-Matveev 
formula for the finite-gap solutions of the KdV equation. 

The formula ( |6.24[ ) derived in |54|, ^| has led to one of the most important pure math- 
ematical applications of the theory of non-linear integrable systems. This is the solution of 
the famous Riemann-Shottky problem. 

According to the Torrelli theorem, the matrix of 6-periods of normalized holomorphic dif- 
ferentials uniquely defines the corresponding algebraic curve. The Riemann-Shottky problem 
is: to describe symmetric matrices with the positive imaginary part which are the matrices 
of 6-periods of normalized holomorphic differentials on algebraic curves. One of the authors 
conjected that the function u(x,y,t) given by fl6.24p is a solution of the KP-equation iff the 



matrix B that defines the theta-function is the matrix of 6-periods of normalized holomorphic 
differentials on an algebraic curve and U, V, W are vectors of 6-periods of corresponding nor- 
malized meromorphic differentials with the only pole at a point of this curve. This conjecture 
was proved in |74jl . 



7 Spectral theory of two-dimensional periodic opera- 
tors. KP hierarchy 

A general algebraic-geometric construction of the finite-gap potentials for the Schrodinger 
operators and for solutions of the KP hierarchy that was presented in the previous section 
has been developed extensively in years. It is applicable for all soliton systems which are 
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equivalent to various types of compatibility conditions for over-determined systems of aux- 
iliary linear problems. In its algebraic form it is in some sense local and is a sort of inverse 
transform: from a set of algebraic-geometrical data to solutions of the integrable non-linear 
partial differential equations 

{algebraic — geometrical data} i — ► {solutions of NLPDE} (7.1) 

In a generic case the space of algebraic-geometrical data is a union for all g of the spaces 

M gtN = {T g ,P a ,k- 1 (Q), ll ,..., lg }, a = l,...,N, (7.2) 

where T g is an algebraic curve of genus g with fixed local coordinates k~ l {Q), k~ 1 (P a ) = 0, 
in neighborhoods of iV punctures P a , and 71, . . . ,7 9 are points of T g in a general position. 
(It is to be mentioned that M g ^ are "universal" data. For the given non-linear integrable 
equation the corresponding subset of data has to be specified.) 

A'posteriory it can be shown that these solutions can be expressed in terms of the corre- 
sponding Riemann theta-functions and are quasi-periodic functions of all variables. Whithin 
this approach it is absolutely impossible to give an answer to the basic question: "How 
many algebraic-geometrical solutions are there? And what is their role in the solution of the 
periodic Cauchy problem for two-dimensional equations of the KP type?" 

The answer of the corresponding question in lower dimensions is as follows. For finite 
dimensional (0+1) systems a typical Lax representation has the form 

d t U(t,\) = [U(t,\),V(t,\)), (7.3) 

where U(t, A) and V(t, A) are matrix functions that are rational (or sometimes elliptic) 
functions of the spectral parameter A. In that case all the general solutions are algebraic- 
geometrical and can be represented in terms of the Riemann theta-functions. 

For spatial one-dimensional evolution equations of the KdV type (l+l)-systems) the 
existence of direct and inverse spectral transform allows one to prove (though it is not 
always the rigorous mathematical statement) that algebraic-geometrical solutions are dense 
in the space of all periodic (in x) solutions. 

It turns out that the situation for two-dimensional integrable equations is much more 
complicated. For one of the real forms of the KP equation that is called the KP-2 equation 
and coincides within ( |6.9| ), the algebraic-geometrical solutions are dense in the space of all 
periodic (in x and y) solutions (|^T|). It seems, that the same statement for the KP-1 equation 
which can be obtained from ( |6.9|) be replacing y — > iy is wrong. One of the most important 
problems in the theory of two-dimensional integrable systems which are still unsolved is: "in 
what sense" the KP-1 equation ahich has the operator representation ( |6.7p and for which a 
wide class of periodic solution was constructed, is an " non-integrable" system. 

The proof of the integrability of the periodic problem for the KP-2 equation is based on 
the spectral Floquet theory of the parabolic operator 

M = d y -d£+u(x,y), (7.4) 
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with periodic potential u(x + h,y) = u(x,y + l 2 ) = u(x,y). We are going now to present 
the most essential points of this theory which was developed in It is the natural 



generalization of the spectral theory of the periodic Sturm-Liouville operator. We would 
like to mention that despite its application to the theory of non-linear equations and related 
topics, the structure of the Riemann surface of Bloch solutions of the corresponding linear 



equation that was found in [41] has been used as a starting point for an abstract definition 



of the Riemann surfaces of the infinite genus ([75]). 

Solutions ip(x,y,Wi,W2) of the non-stationary Shrodinger equation 

(ad y - dl + u(x, y))ip{x, y, wi, w 2 ) = (7.5) 

with a periodic potential u(x,y) = u(x + a\,y) = u(x,y + a 2 ) are called Bloch solutions, if 
they are eigenfunctions of the monodromy operators, i.e. 

ip(x + a x ,y, wi, w 2 ) = Wiip(x, y, wi, w 2 ), (7.6) 

ip(x, y + a 2 , wi, w 2 ) = w 2 tp(x, y, wi, w 2 ). (7.7) 

The Bloch functions will always be assumed to be normalized so that ^(0, 0, w\, w 2 ) = 1. 
The set of pairs Q = (wi,w 2 ), for which there exists such a solution is called the Floquet set 
and will be denoted by T. The multivalued functions p(Q) and E(Q) such that 

Wl = e ipai , w 2 = e iEa2 

are called quasi-momentum and quasi-energy, respectively. 

The gauge transformation ip — > e^^ip, where d y h{y) is a periodic function, transfers the 
solutions of ( [7.5p into solutions of the same equation but with another potential u = u — od y h. 
Consequently, the spectral sets corresponding to the potentials u and u are isomorphic. 
Therefore, in what follows we restrict ourselves to the case of periodic potentials such that 
Jq 1 u(x, y)dx = 0. 

To begin with let us consider as a basic example the "free" operator Mo = crd y — dl with 
zero potential u(x,y) = 0. The Floquet set of this operator is parameterized by the points 
of the complex plane of the variable k = e tkai , w 2 = e~ u k 0,2 , and the Bloch solutions 
have the form ip(x,y,k) = e tkx ~ u k y . The functions tp + (x,y,k) = e~ lkx + a k v are Bloch 
solutions of the formal adjoint operator (crd y + dl)ip + = 0. 

An image of the map k e C i — > (w^w®) is the Floquet set for the free operator M . It 
is the Riemann surface with self-intersections. The self-intersections correspond to the pairs 
k 7^ k! such that wfik) = wf(k'), i — 1, 2. The latter conditions imply the equations 

k — k — , k - (k'f = , (7.8) 

d a 2 

where iV and M are integers. Hence, all the resonant points have the form 

nN aiMa\ , . . 

k = k NM = — N^0, k — k- N) - M . 7.9 

ai JMa 2 
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The basic idea of the construction of the Riemann surface of Bloch solutions of the equation 
( [7.5[) that was proposed in |4l| is to consider (|7.5| ) as a perturbation of the free operator, 



assuming that the potential u(x,y) is formally small. 

For any k ^ k N ^ M it is easy to construct a formal Bloch solution. It turns out that the 
corresponding formal series converges and defines a holomorphic function of k for \ko\ > M 
big enough and lies outside small neighborhoods of the resonant points. Moreover, it can be 
shown that this function can be extended on the Riemann surface that can be thought as a 
surface obtained from the complex plane by some kind of surgery that creates gaps in places 
of the resonant points. 

More precisely, if u(x,y) is a smooth real potential that has analytical continuation in 
some neighborhood of the real values of x and y, then the corresponding Riemann surface 
of Bloch-Floquet solutions can be described in the following way. 

Let us fix some finite or infinite subset S of integer pairs (N > 0, M) . The set of pairs 
of complex numbers ix = {p s ,i,p s ,2} where s G S would be called "admissible", if 

Re p s ,i = — , \p s ,i ~ K\ = oQ^r 1 ), i = 1, 2, (7.10) 

and the intervals \p s ,i,p s ,2] do not intersect. (Here k s , s = (N,M) are resonant points.) 

Let us define the Riemann surface r(7r) for any admissible set tt. It is obtained from the 
complex plane of the variable k by cutting it along the intervals \p s ,i,Ps,2] and [— p s ,i, — p s ,2] 
and by sewing after that the left side of the first cut with the right side of the second cut 
and vice versa. ( After this surgery for each cut [p s ,i,p s ,2\ corresponds to a nontrivial cycle 
a s on r(7r).) 

For any real periodic potential u(x, y) which can be analytically extended into some 
neighborhood of the real values x, y, the Bloch solutions of the equation ( [7.5|) are parame- 
terized by points Q of the Riemann surface T(tt) corresponding to some admissible set 7r. 
The function ip(x, y, Q) which is normalized by the condition ip(0, 0, Q) — 1 is meromorphic 
on T and has a simple pole 7 S on each cycle a s . If the admissible set it contains only a 
finite number of pairs, then T(7r) has finite genus and is compactified by only one point P\ 
(k = oo), in the neighborhood of which the Bloch function ip has the form ( |6.2| ). 

The potentials u for which r(7r) has finite genus are called finite-gap. They coincide with 
the algebraic-geometrical potentials. The direct spectral transform for periodic operators 
( [7.4| ) allows us to prove that like in the one-dimensional case the finite-gap potentials are 
dense in the space of all periodic smooth functions in two variables (||41||). 



8 Spectral theory of two-dimensional Schrodinger op- 
erator for fixed energy level and two-dimensional 
Toda lattice 

In this section we discuss a spectral theory of two-dimensional periodic Schodinger operator. 
Unlike the one- dimensional case, spectral data for 2D linear operator are over-determined 
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and therefore for generic operators there are no nontrivial isospectral flows. As it was noted 
in 1 77[1 , deformations that preserve spectral data for one fixed energy level do exist. An 



analog of Lax representation for such system has the form 

H t = [A, H] + BH, (8.1) 

where H, A, B are two dimensional operators with coefficients depending on x, y, t. Equation 
( |3.1D is equivalent to the condition that operators H and (d t — A) commute on the space of 
solutions of the equation Hip = 0. Therefore, ( |8.1| ) describes deformations preserving all the 
spectral data associated with the zero energy level of the operator H. 



It should be mentioned that until the moment when equations (8.1) were proposed in 
the framework of the soliton theory the spectral problem associated with one energy level of 
two-dimensional periodic operators had never been considered. 

For the first time an inverse algebraic-geometric spectral problem for two-dimensional 
Schrodinger operator in magnetic field 

H = (id x - A x {x, y)f + (id y - A y (x, y) 2 + u(x, y) (8.2) 



was formulated and solved in 36 



Consider the Bloch solutions of the equation Hip = eip, which by definition are solutions 
that at the same time are eigenfuctions for the shifts operators: 

iP(x + 7i, y) = e ip * Tl iP{x, y), if>(x, y + T 2 )= e ip ^(x, y). (8.3) 

Here T 1; T 2 are periods of the operator H, i.e. periods of the potential u(x, y) and periods of 
the magnetic field B(x, y), which is defined by the vector potential (A x , A y ) by the formula 
B = dyA x d x Ay. 

Multivalued quantities p x ,p y are components of the two-dimensional quasimomentum. 
For fixed values of p x ,p y a spectrum of the operator H restricted on the space of functions 
satisfying ( |8.3| ) is discrete and defines different branches of dispersion relations €j{p x ,p y ),j = 
1, . . . Level lines €j(p x ,p y ) = e in the space of variables p x ,p y define the, so-called, Fermi 
curves. Of course, in the solid state physics all the considerations were restricted by real 
values of quasi-momentum. 



In [ 36 1 it was suggested to consider operators for which a complex Fermi curve does 
exist and is the Riemann surface of finite genus for some energy level eo- Moreover, it was 
assumed that this curve is compactified by two infinity points P± in neighborhoods of which 
the corresponding Bloch solutions have the form 

4 = e^^ ) (£g{x,y)k?y (8.4) 

where A;^ 1 ) are local coordinates in the neighborhoods of the punctures P±. It was also 
assumed that outside the punctures the functions ip(x,y,P) considered as a function of the 
variable P (which is a point of the complex Fermi curve T) is meromorphic function with g 
poles independent of the variables x and y. 
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We present here a solution to this inverse spectral problem in a more general form which 
is necessary for construction of exact solutions to the two-dimensional Toda lattice which 
has deep connections with the theory of 2D Schrodinger operators. After that we return 
back to the spectral problems. 

Let T be a smooth algebraic curve of genus g with fixed local coordinates w±(P) in 
neighborhoods of the points P ± , w±(P ± ) = 0. Then for each set of g points 71, . . . , 7 5 in 
general position there exists a unique function ip n (T : P), T = {tf, % — 1, . . . , 00, } such that: 

1°. The function ip n of the variable P E T is meromorphic outside the punctures and has 
at most simple poles at the points 7 S (if all of them are distinct); 

2°. In a neighborhood of the point P^ 1 it has the form 



1>- 



n (T,P) = wf 1 {^if{n,T)w s ^ exp (X>±^ , w ± = w ± (P), (8.5) 

ti(x,T) = 5 ar (8.6) 

The proof of this statement, as well as the explicit formula for ip n in terms of Riemann 
theta-functions, is almost identical to the way of solution of the inverse problem for finite- 
gap Schrodinger operators discussed in the previous sections. 

Let dVt^ be a unique meromorphic differential holomorphic on V outside the point P ± , 
which has the form 

dVtf ] = d(wZ 3 + 0(w±)), (8.7) 



3 

near the point P ± , and normalized by the conditions § dVt^ = 0. It defines a vector 
with coordinates 

1 

Further, let us define the normalized differential dVL^\ which is holomorphic outside the 
points P^ where it has simple poles with residues ±1, respectively. From Riemann's bilinear 
relations it follows that the vector of 6-periods of this differential equals 2niU^ , where 

C/(°) = A(P-) - A(P + ) . (8.9) 

As in the previous case, one can directly check that the function ip n {T, P) given by the 
formula: 

4(TP) = 9{A{P) + U( ° )n + — U ^ tf + Z )°^ P+ ) + g) e (^(m£Wp)) (8 10) 
' 6{A{P) + Z)d{A{P+) + U^n + Y J U ( i ±) tf + Z) 

a (±) (p) = fdn^ , (8.11) 

is well-defined and has all the properties of the Baker- Akhiezer function. 

Note, that for n = and tf — x ± iy, tf — 0, % > 1, the analytical properties of this 
function coincide with the properties that were described above as analytical properties of 
the Bloch solutions for finite-gap two-dimensional Schrodinger operators. 
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^From the uniqueness of the Baker- Akhiezer functions il) n (T, P) it follows that they satisfy 
the linear equations 

d 

d + i) n = l[) n+1 + V n l/>n, = Cn^n-l, 9± = 7^, (8.12) 

where 

v n = d + <p n (T), Cn = e^-^ T \ e v " = £o(T), (8.13) 

and £,o(T) is a leading term of the expansion of ip n at the puncture P~ . From (|8.10| ) we get 
the formula 

i 0{U {Q) {n + l) + Y.U { i ±) tt + Z Q ) 7 ^Aiv+\ / 8 ,/i 

ip n = log 7-t h const, Zn = Z + A(P^). 8.14 

for algebraic-geometric solutions of 2D Toda lattice which was obtained in |76 
Note that Q^.12 ) imply that ?/>o satisfies the equation 

d+d-ipo + v d-i/j + (ci + d-v )i(>o, (8.15) 



which is gauge equivalent to ( ft.2|) and therefore, we do get a solution of the inverse problem 
that was introduced above. 

The next important step was done in |38], |39| where algebraic-geometric spectral data 
corresponding to potential 2d Schrodinger operators (i.e. operators with zero magnetic field) 
were found. 

Let T be a smooth genus g algebraic curve with fixed local coordinates k^}{Q), 
k^(P±) = 0, in neighborhoods of two punctures P±. Let us assume that there exists a 
holomorphic involution of the curve a : V i — > V such that P± are its only fixed points, 
i.e. cr(P±) = P±. The local parameters are to be "odd", i.e. k±(cr(Q)) = —k±(Q). The 
factor-curve will be denoted by r . The projection 

7T : T i — > T = Y/a (8.16) 

represents T as a two-sheet covering of To with the two branch points P±. In this realization 
the involution a is a permutation of the sheets. As there are only two branching points 
g = 2g , where go is genus of T . Let us consider a meromorphic differential dQ(Q) of the 
third kind on r with residues =pl at the points P±. The differential dVL has g zeros that will 
be denoted by i = 1, . . . , 2g = g. Let us choose for each i a point ji on T such that 

7r(7i) =li, i = l,...,g. (8.17) 



In |38], |39| it was shown that the Baker- Akhiezer function corresponding to algebraic-geomet- 
ric data which have been just defined satisfies the equation 

(d + d_ + u(x,y))^(x,y,Q) = 0, (8.18) 

where 

u = -d-z+ = -a+er, (8.i9) 
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and ^ = ^ (x + , X-) are the first coefficients in the expansion ( |S.4|) . 

It should emphasized that though general formula ( |3.10|) in terms of the Riemann 9- 
functions is valid for the Baker- Akhiezer functions corresponding to the potential Shrodinger 



operator H, in [38, 39] it was found another more effective representation in terms of the, 



so-called, Prim theta-f unction. 

A space of holomorphic differentials on T splits into two go-dimensional subspaces of 
even and odd (with respect to the involution a) differentials. A matrix of 6-periods of odd 



differentials defines the function 9p r (z) by the same formula ( 6.14 ). Then 



V 6 Pr {A° d {Q) - Z)6 Pr (U+x+ + U-x- -Z) ' K ' 

Here p (Q) are Abelian integrals of the second kind normalized differential dp* 1 on V that 
have poles of the second order at the points P±, respectively; vectors 2irU ± are the vectors 
of 6-periods of these differentials. 

As it was mentioned above the inverse algebraic-geometric spectral problem on one energy 
level for 2D Schrodinger operator was posed and solved at the time when no direct spectral 
theory was known. This theory was developed much later in jy]], where it was shown that 
the Bloch solutions for ( |8.18| ) with analytical periodic potential are parameterized by points 
of infinite genus Riemann surface. It was proved that if this surface has finite genus, then 
the Bloch functions have all the analytical properties suggested in the inverse problem and 
therefore are just the Baker-Akhiezer-functions. Moreover, it was proved that algebraic- 
geometric (finite-gap) potentials are dense in the space of all periodic potentials. 



9 Spectral theory of operators with elliptic coefficients 

In this section we are going to discuss a specific spectral problem for operators with elliptic 
coefficients, i.e. with the coefficients that are meromorphic functions of a variable x and 
have two periods 2u, 2lu', Im (uj'/uj) > 0. 

Since Hermit time, it has been known that one- dimensional Schrodinger operator with 
potential of the form n(n + l)p(x) (where p{x) = p(x\u,u') is Weierstrass p-function corre- 
sponding to elliptic curve with periods 2u, 2u', and n is an integer) has only n gaps in the 
spectrum. These Lame potentials had been the only known examples with the finite-gap 
property before the finite-gap theory was constructed in the framework of the soliton theory 
(see above). As we already have shown a generic algebraic-geometric potential can be ex- 
pressed in terms of higher genus Riemann ^-function. Sometimes higher genus formula can 
be reduced to exact expression in terms of elliptic function. 

The first example of this type which is different from the Lame potentials was proposed 
in [25]. Later a theory of elliptic finite-gap potentials attracted the special interest due to 



remarkable observation made in [78] on their connection with the elliptic Calogero-Moser 
model. The most recent burst of interest is due to the unexpected connections of these 



systems to Seiberg-Witten solution of N = 2 supersymmetric gauge theories [14, 15]. It 
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turns out that the low energy effective theory for SU(N) model with matter in the adjoint 
representation (identified first in |l| with SU(N) Hitchin system) is isomorphic to the elliptic 
CM system. Using this connection quantum order parameters were found in |79f. 



The elliptic Calogero-Moser (CM) system ([80, pl]| ) is a system of N identical particles 
on a line interacting with each other via the potential V(x) = p(x). Its equations of motion 
have the form 

The CM system is a completely integrable Hamiltonian system, i.e. it has N independent 
integrals Hk in involution ( ||82|| , ||S"3"|). The second integral if 2 is the Hamiltonian of ( |9.1Q . 



In [78] it was shown that the elliptic solutions of the KdV equations have the form 

N 

u(x,t) = 2j2p(x-Xi(t)) (9.2) 



i=l 



and the poles Xi(t) of the solutions satisfy the constraint J2j^i p'( x i ~ x j) = 0, which is the 
locus of the stationary points of the CM system. Moreover, it turns out that the dependence 
of the poles with respect to t coincides with the Hamiltonian flow corresponding to the 
third integral H s of the system. In |8~§, B3] it was found that this connection becomes 



an isomorphism in the case of the elliptic solutions of the Kadomtsev-Petviashvilii equation. 
Moreover, in |6(| it was revealed that the connection of the CM systems with the KP equation 



is in some sense secondary and is a corollary of more fundamental connections with spectral 
theory of linear operators with elliptic potentials. The corresponding approach has been 
developed extensively in |^7|, ^ |55| . 

Let £ be a linear differential or difference operator in two variables x, t with coefficients 
which are scalar or matrix elliptic functions of the variable x. We do not assume any special 
dependence of the coefficients with respect to the second variable. Then it is natural to 
introduce a notion of double-Bloch solutions of the equation 

Cip = 0. (9.3) 

We call a meromorphic vector-function f(x), which satisfies the following monodromy prop- 
erties: 

f{x + 2u a ) = B a f(x), a = 1,2, (9.4) 

a double-Bloch function. The complex numbers B a are called Block multipliers. (In other 
words, / is a meromorphic section of a vector bundle over the elliptic curve.) 

In the most general form a problem that we are going to address is to classify and to 
construct all the operators £ such that equation (|9.3| ) has sufficiently enough double-Bloch 
solutions. 

It turns out that existence of the double-Bloch solutions is so restrictive that only in 
exceptional cases such solutions do exist. A simple and general explanation of that is due to 
the Riemann-Roch theorem. Let D be a set of points Xj, % = 1, . . . , m, on the elliptic curve 
To with multiplicities di and let V = V(D; Bi, B 2 ) be a linear space of the double-Bloch 
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functions with the Bloch multipliers B a that have poles at Xi of order less or equal to di and 
holomorphic outside D. Then the dimension of D is equal to: 

dim D = deg D — dj. 

i 

Now let Xi depend on the variable t. Then for / e D{t) the function Cf is a double- 
Bloch function with the same Bloch multipliers but in general with higher orders of poles 
because taking derivatives and multiplication by the elliptic coefficients increase these orders. 
Therefore, the operator £ defines a linear operator 

C\ D : V(D(t);B 1 ,B 2 ) .— > V(D'(t); B x , B 2 ), N 1 = degD' > N = degD, 

and ( |9.3| ) is always equivalent to an over- determined linear system of N' equations for iV 
unknown variables which are the coefficients q = Ci(t) of expansion of ^ G V(t) with respect 
to a basis of functions fi(t) e V(t). With some exaggeration one may say that in the 
soliton theory the representation of a system in the form of the compatibility condition of 
an over-determined system of the linear problems is considered equivalent to integrability. 

In all of the basic examples N' = 2N and the over-determined system of equations has 
the form 

LC = kC, d t C = MC, (9.5) 

where L and M are N x N matrix functions depending on a point z of the elliptic curve as 
on a parameter. A compatibility condition of ( |9.5|) has the standard Lax form d t L = [M, L], 
and is equivalent to a finite-dimensional integrable system. 

The basis in the space of the double-Bloch functions can be written in terms of the 
fundamental function z) defined by the formula 

= (9.6) 

Note, that <&(x, z) is a solution of the Lame equation: 



^ - 2p(x) ) $(x, z) = p(z)$(x, z). (9.7) 



^From the monodromy properties it follows that $ considered as a function of z is double- 
periodic: 

$(x,z + 2u a ) = $(x,z), 

though it is not elliptic in the classical sense due to an essential singularity at z = for 
x ^ 0. 

As a function of x the function <&(x, z) is double-Bloch function, i.e. 

$(x + 2u a , z) = T a (z)$(x, z), T a (z) = exp (2u a ((z) - 2((uj a )z) . 

In the fundamental domain of the lattice defined by 2u a the function z) has a unique 
pole at the point x = 0: 

$(x,z) = x' 1 + 0(x). (9.8) 
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The gauge transformation f(x) i — > f{x) = f(x)e ax , where a is an arbitrary constant does 
not change poles of any function and transforms a double Bloch-function into another double- 
Bloch function. If B a are Bloch multipliers for / then Bloch multipliers for / are equal to 

S 1 = S 1 e aaUl , B 2 = B 2 e 2aW2 . (9.9) 

The two pairs of Bloch multipliers that are connected with each other through the relation 
( p.9[) for some a are called equivalent. Note that for all equivalent pairs of Bloch multipliers 
the product B^B^ 1 is a constant depending on the equivalence class, only. 

^From ( |9.8| ) it follows that a double-Bloch function f(x) with simple poles Xi in the 
fundamental domain and with Bloch multipliers B a (such that at least one of them is not 
equal to 1) may be represented in the form: 

N 

f(x) = '£c i <S>{x-x i ,z)e kx , (9.10) 

where c, is a residue of / at X{ and Jv £1X6 parameters related by B a = T a (z)e 2uJak . (Any 
pair of Bloch multipliers may be represented in this form with an appropriate choice of the 
parameters z and k.) 

Let us consider as an example the equation 

&/) = (d t -c% + u(x,t))ip = 0, (9.11) 



where u(x,t) is an elliptic function. Then as shown in equation ( |9.11| ) has N linear 



independent double-Bloch solutions with equivalent Bloch multipliers and N simple poles at 
points Xi(t) if and only if u(x,t) has the form 

N 

u{x,t) = 2j2p{x-Xi(t)) (9.12) 
i=i 

and Xi(t) satisfy the equations of motion of the elliptic CM system ( p.l| ). 

The assumption that there exist N linear independent double-Bloch solutions with equiv- 
alent Bloch multipliers implies that they can be written in the form 

N 

i> = J2 ~ Xi®, z)e kx+k2 \ (9.13) 

i=i 

with the same z but different values of the parameter k. 



Let us substitute (|9.13|) into ( |9.11| ). Then ( |9.11| ) is satisfied if and only if we get a function 



holomorphic in the fundamental domain. First of all, we conclude that u has poles at x« 
only. The vanishing of the triple poles (x — Xj)~ 3 implies that u(x,t) has the form ( |9.12| ). 



The vanishing of the double poles (x — Xj) 2 gives the equalities that can be written as a 
matrix equation for the vector C = (q): 

(L(t, z) - kI)C = , (9.14) 
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where I is the unit matrix and the Lax matrix L(t, z) is defined as follows: 

L v (t,z) = -±6 ij i i -(l-6 ij mx i -x j ,z). (9.15) 
Finally, the vanishing of the simple poles gives the equations 

(d t -M(t,z))C = 0, (9.16) 

where 

Ma = (p{z) - 2 P( x i ~ x j)\ ~ 2(1 - 6y)&{xi ~ x jt z). (9.17) 



The existence of N linear independent solutions for (|9.11|) with equivalent Bloch multipliers 



implies that ( |9.14| ) and fl9.16| ) have N independent solutions corresponding to different values 
of k. Hence, as a compatibility condition we get the Lax equation L = [M, L] which is 
equivalent to ( |9.1[ ). Note that the last system does not depend on z. Therefore, if fl9.14| ) 



and (|9.16|) are compatible for some z, then they are compatible for all z. As a result we 
conclude that if ( |9.11| ) has N linear independent double-Bloch solutions with equivalent 
Bloch multipliers then it has infinitely many of them. All the double-Bloch solutions are 
parameterized by points of an algebraic curve F defined by the characteristic equation 



N 



R(k, z) = det(kl - L{z)) = k N + J2 niz)^-* = 0. (9.18) 



Equation ( p .18 ) can be seen as a dispersion relation between two Bloch multipliers and 



defines T as iV-sheet cover of IV 



As was shown in [66] expansion of the characteristic equation ( 9.14 ) at z = has the 
form: 

N 

R(k,z) = H(k + is i z- 1 + hi + 0(z)), v x = 1 - JV, u { = 1, i > 1. (9.19) 

i=l 

We call the sheet of V at z = corresponding to the branch k = z~ x (N — 1) + 0(1), an 
upper sheet and mark the point P\ on this sheet among the preimages of z = 0. From ( p.!9| ) 
it follows that in general position when the curve T is smooth, its genus equals N. 

Further consideration of analytical properties of the function tp given by fl9.13|) where 
Cj are components of the eigenvector to the matrix L shows that this function is just the 
Baker- Akhiezer function introduced in Section 6. Combined with expression ( |6.22| ) for if) in 
terms of ^-function this result directly leads to the main statement of |6l| : 

the coordinates of the particles Xi(t) are roots of the equation 

9(Ux + Vt + Z Q ) = 0, (9.20) 

where = 8(£\B) is the Riemann theta-function corresponding to the matrix of 6-periods 
of holomorphic differentials on T; the vectors U and V are the vectors of 6-periods of nor- 
malized meromorphic differentials on T, with poles of order 2 and 3 at the point Pi. 
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Among other examples of integrable systems that can be generated in the similar way 
are Ruijesenaars- Schneider system |86 



.1, j 



XiX s (V(xi - x s ) - V(x s - Xi)), V(x) = C(x) - ((x + 77), 



(9.21) 



and nested Bethe ansatz equations |87| 



a(x™ - x™ +1 )a(xt —T} — x")a(a;" - x™ 1 + rj) 



-1 



As shown in 



and 



]J a(x? - rif L )<r{x* + r]~ x])a(x? - x] +1 - rj) 

they are generated by spectral problems for equations 
Cip = d t ip(x, t) — ip(x + r),t) — v(x, t)ip(x, t) = 0, 



(9.22) 

(9.23) 
(9.24) 



-0(x, m + 1) = vp(x + rf) + v(x, m)^{x } m), 
respectively. (Here rj is a complex number and v(x,t) is an elliptic function.) 

Strange as it is, the inverse spectral problem which is discussed here is simplier for two- 
dimensional operators then for one-dimensional stationary operators. For example a family 



of spectral curves corresponding to operators (|9.11|) that have double-Bloch solutions can be 
described explicitly. A nice formula was found in |79| : 



(9.25) 



where H(k) is a polynomial. Note that ( p.25 ) may be written as: 



/(M 



N 1 



dk 



H(k). 



The coefficients of the polynomial H (k) are free parameters of the spectral curve of the CM 
system. 

The spectral curves corresponding to the Schrodinger operator with the same property 
are special case of the curves (|9.25|) but their explicit description is unknown. In particular, 
exact formula for branching points for Lame potentials is unknown. As it was mentioned 
above the first example of elliptic finite-gap potentials different from the Lame potentials 
was found in ||25|| . A wide class of such potentials was found in ||70 . 

In |88| it was noted that the problem of classification of Schrodinger operators with elliptic 
potentials that have two double-Bloch solutions for almost all energy levels was posed by 
Picard though had not been solved until very recently. In [R3 using Floquet spectral theory 



for the Schrodinger operator it was proved that all such potentials are finite-gap. This result 
is an essential step in the Picard problem though its complete and effective solution is still 
an open problem. 
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